PROCEEDINGS  of  the  THIRD 
BERKELEY  SYMPOSIUM  ON 
MATHEMATICAL  STATISTICS 
AND  PROBABILITY 

Held  at  the  Statistical  Laboratory 
University  of  California 
December  26-37,  1954 
C'July  and  August \  7955f) 


VOLUME  II 

Contributions  to  Probability  Theory 
Edited  by  JERZY  NEYMAN 


For  further  effective  support  of  the  Symposium  thanks  must  be  given  the 
National  Science  Foundation,  the  United  States  Air  Force  Research  and  Development 
Command,  the  United  States  Army  Office  of  Ordnance  Research,  and  the  United  States 
Navy  Office  of  Naval  Research. 


UNIVERSITY  OF  CALIFORNIAJPRESS 

BERKELEY  AND  LOS  ANGELES 


1956 


THE  ZEROS  OF  A  RANDOM 
POLYNOMIAL 


J.  M.  HAMMERSLEY 

OXFORD  UNIVERSITY 


1.  Introduction 


If  the  coefficients  of  an  algebraic  equation  are  subject  to  random  error,  the  roots  of 
this  equation  will  also  be  subject  to  random  error;  and  it  is  natural  to  enquire  how  the 
latter  errors  depend  upon  the  former.  This  is  clearly  a  question  of  some  practical  impor¬ 
tance.  It  arises  not  only  when  the  coefficients  result  from  experimental  data,  but  also, 
for  example,  when  the  coefficients  are  rounded  off  to  some  specified  number  of  decimal 
places  before  commencing  a  numerical  solution.  Yet  it  is  a  question  which  has  so  far  re¬ 
ceived  rather  scant  attention,  apart  from  the  treatment  of  three  special  instances. 

In  the  first  of  these  special  instances  the  equation  is  of  a  particular  type,  namely  the 
characteristic  equation  of  a  variance-covariance  matrix  pencil  whose  elements  are  real 
and  distributed  in  Wishart’s  form.  Under  these  circumstances  the  roots  are  all  real,  and 
their  joint  sampling  distribution  is  well  known.  Amongst  the  several  textbooks,  which 
discuss  this  question,  the  reader  may  consult  Wilks  [19]. 

In  the  second  special  instance  the  equation  is  linear  and  the  coefficients  are  real  and 
distributed  normally  (though  not  necessarily  independently).  The  distribution  of  the 
root  of  this  equation  is  therefore  that  of  the  quotient  of  two  real  correlated  normal  vari¬ 
ates.  Geary  [9]  gives  the  required  result.  This  special  instance  arises  in  bio-assay  work 
under  the  name  of  Fieller’s  theorem  (see  pp.  27-29  in  Finney  [8]).  The  interpretation 
of  this  theorem  is,  however,  open  to  question,  and  I  discuss  this  matter  further  in  sec¬ 
tion  9. 

In  the  third  special  instance  the  equation  is  of  general  degree  and  its  coefficients  are 
real  and  distributed  independently  and  symmetrically  about  zero  either 

(i)  normally,  or 

(ii)  rectangularly,  or 

(iii)  discretely  into  the  pair  of  classes  ± 1. 

Littlewood  and  Offord  [14],  [15]  enquired  how  many  real  roots  on  the  average  such  a 
random  equation  might  be  expected  to  have,  and  they  gave  asymptotic  approximations 
for  the  result  valid  for  equations  of  large  degree.  Kac  [12],  [13]  improved  their  results 
by  showing  that  the  average  number  of  real  roots  of  an  equation  of  degree  n  —  1  is 
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when  (i)  applies;  and  for  both  (i)  and  (ii)  Nn  ~  ( 2/x )  log  n  as  n— >  °o.  For  the 
latest  results,  see  Erdos  and  Offord  [20].  Littlewood  and  Offord  noticed  that  the 
three  cases  (i),  (ii),  and  (iii)  led  to  results  which  were  (as  far  as  they  could  tell 
from  their  approximations)  all  of  the  same  order  of  magnitude  for  large  n.  This  is 
not  surprising  in  view  of  the  result  obtained  in  section  6  that  the  average  number 
of  real  roots  depends  only  on  the  distribution  of  the  real  and  imaginary  parts  of  the 
polynomial  and  its  derivatives  (regarding  its  argument  as  fixed);  for  by  the  central 
limit  theorem  these  four  parts  will  (under  wide  conditions)  be  asymptotically  normal. 
However,  as  Kac  [13]  discovered,  it  is  not  easy  to  put  this  reasoning  into  rigorous  terms; 
and  we  shall  not  attempt  it  here. 

Some  quotations  from  the  papers  of  Littlewood  and  Offord  and  of  Kac  indicate  certain 
differences  of  approach  in  the  present  article.  From  [15]:  “In  ‘selecting’  a  random  equa¬ 
tion  of  a  particular  given  degree  n  (for  example,  n  =  10)  it  is  natural  to  treat  each  coeffi¬ 
cient  on  the  same  basis.”  This  is  a  condensation  of  a  longer  passage  in  [14],  where  two 
players  are  betting  on  the  number  of  real  roots  of  an  equation  selected  by  a  referee,  who 
“would  presumably  treat  each  coefficient  on  the  same  basis.”  Similarly  from  [13]:  “.  .  . 
let  [the  coefficients]  be  independent  random  variables  each  having  a(u)  as  its  distribu¬ 
tion  function.”  It  is,  of  course,  reasonable  to  treat  the  coefficients  as  independent  and 
having  the  same  distribution  function  when  the  problem  is  posed  as  one  of  pure  mathe¬ 
matics;  and,  even  when  entering  upon  a  wider  sphere  than  that  of  pure  mathematics, 
it  is  doubtless  wise  to  regard  the  simple  case  of  similarly  distributed  coefficients  as  a  start¬ 
ing  point  in  one’s  investigations.  But  in  a  problem  of  applied  mathematics  it  is  far  from 
likely  that  the  coefficients  will  be  independent  or  have  similar  distributions.  My  own  in¬ 
terest  in  the  problem  sprang  from  a  question  put  to  me  concerning  the  precision  of  an 
estimate  of  the  growth  rate  of  an  insect  population,  which  growth  rate  is  the  root  of  an 
algebraic  equation  whose  coefficients  are  determined  from  experimental  data  in  a  manner 
that  certainly  vitiates  the  assumption  of  independence.  So  in  the  present  article  I  allow 
the  coefficients  to  be  distributed  in  a  quite  general  manner.  Although  I  have  been  able 
to  determine  a  formal  expression  for  the  distribution  of  the  roots  in  this  general  case,  I 
am  still  very  far  from  having  solved  the  practical  problem;  for  my  result  involves  the 
parameters  of  the  distribution  function  of  the  coefficients,  and  except  in  a  trivial  special 
case  I  cannot  at  present  see  how  to  “Studentize”  the  result  in  such  a  way  that  these 
parameters  may  be  replaced  by  their  estimates.  Concerning  this  problem  of  “Studentiza- 
tion,”  the  reader  may  care  to  consider  the  papers  by  Creasy  [4]  and  Fieller  [7]. 

Another  result  of  attacking  the  problem  from  the  point  of  view  of  generally  distributed 
coefficients  was  to  throw  light  upon  certain  essential  points  of  difficulty  in  the  problem. 
In  [13]  Kac  writes:  “Upon  closer  examination  it  turns  out  that  the  proof  I  had  in  mind 
[in  [12]]  ...  is  inapplicable  to  the  [case  when  the  coefficients  have  a  discrete  distribu¬ 
tion].  .  .  .  This  situation  tends  to  emphasize  the  particular  interest  of  the  discrete  case, 
which  surprisingly  enough  turns  out  to  be  the  most  difficult.”  Something  more  than  this 
is  in  fact  surprising:  namely,  that  the  simplest  case  arises  when  the  coefficients  are  al¬ 
lowed  to  assume  complex  values,  the  joint  distribution  of  all  real  and  imaginary  parts 
having  a  continuous  frequency  function  with  finite  moments.  When  any  discontinuities 
are  introduced  into  the  distribution  function  of  the  coefficients,  the  corresponding  distri¬ 
bution  of  the  roots  is  in  danger  of  acquiring  discontinuities  as  well.  Thus  when  all  real 
and  imaginary  parts  are  distributed  normally  with  strictly  positive  variances  (that  is 
to  say  in  nondegenerate  multivariate  Gaussian  form)  the  distribution  of  the  roots  pos- 
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sesses  a  frequency  function  over  the  complex  plane;  but  when  the  coefficients  are  made 
real  (so  that  the  distribution  function  of  the  imaginary  part  contains  a  step  function)  the 
distribution  function  of  the  roots,  instead  of  being  a  well-behaved  surface  distribution  in 
the  complex  plane,  partly  gathers  itself  up  into  a  line  distribution  on  the  real  axis,  so 
that  a  frequency  function  no  longer  exists  over  the  whole  complex  plane.  The  results  I 
have  just  quoted  upon  normally  distributed  coefficients  will  be  discussed  in  sections  8 
and  9.  Also  in  the  case  of  normally  distributed  coefficients,  I  shall  fulfill  the  requirement 
in  [12]:  .  .  the  problem  of  the  exact  determination  of  the  average  distribution  of  real 

roots  on  the  real  axis  will,  of  course,  depend  on  a  more  delicate  treatment.” 

In  all  these  three  special  instances,  cited  above,  the  preoccupation  is  with  real  roots. 
As  soon  as  we  attempt  to  tackle  the  more  general  enquiry,  which  embraces  complex 
roots  as  well,  we  come  up  against  a  difficulty.  In  order  to  give  the  joint  distribution  of  a 
set  of  roots  it  is  necessary  to  keep  track  of  each  individual  root.  When  the  roots  are  real 
we  may  do  this  (as  is  done  in  the  first  special  instance)  by  arranging  them  in  order  of 
magnitude.  However,  while  it  is  possible  to  arrange  complex  numbers  in  an  order  of 
magnitude  by  means  of  appropriate  conventions,  it  seems  difficult  to  arrange  them  in  an 
order  which  is  continuous  with  respect  to  the  geometry  of  the  complex  plane.  It  is  con¬ 
ceivable  that,  by  some  device  or  other,  one  might  manage  to  label  individual  complex 
roots  in  such  a  manner  that  their  resulting  joint  distribution  possesses  a  coherent  and 
useful  interpretation;  but  I  have  not  succeeded  in  doing  this,  and  I  am  indeed  doubtful 
whether  it  is  possible.  Alternatively,  one  may  abandon  the  conventional  idea  of  a  joint 
distribution.  In  this  article,  I  have  replaced  it  by  a  somewhat  similar  concept,  which  I 
shall  call  a  condensed  distribution.  We  shall  see  in  due  course  that  this  concept  leads  to 
a  particularly  simple  solution  of  the  problem. 

Another  branch  of  investigation,  connected  with  the  present  problem,  is  followed  by 
Erdos  and  Turan  [5].  They  consider  the  uniformity  of  the  distribution  of  the  roots;  they 

n 

prove  that,  if  N(a,  /3)  denotes  the  number  of  zeros  of  g(z)  =  CjZi  =  0  which  lie  in 

;=0 

the  sector  0  ^  a  <  arg  z  <  /3  ^  2r,  then 


(1.3) 
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2.  Notation 

Capital  German  letters  will  denote  finite-dimensional  Euclidean  spaces.  Small  Ger¬ 
man  letters  and  small  boldface  letters  will  respectively  denote  subsets  and  points  (col¬ 
umn  or  row  vectors)  of  these  spaces.  Capital  boldface  letters  will  denote  matrices.  Italic 
letters  will  denote  scalars  (real  or  complex  according  to  the  context).  It  will  be  con¬ 
venient  sometimes  to  regard  a  complex  scalar  as  a  two-dimensional  vector;  and  we  utilize 
the  foregoing  conventions  by  writing  z  =  x  +  iy  or  z  ={x,  y]  according  to  require¬ 
ments.  A  dagger  attached  to  a  boldface  letter  will  denote  the  transposed  matrix  or  vec¬ 
tor.  A  prime  attached  to  a  German  letter  will  denote  the  complementary  set  (with  re¬ 
spect  to  the  Euclidean  space  containing  the  set  in  question),  while  a  prime  attached  to 
an  italic  or  boldface  letter  will  denote  a  derivative.  |  x  |  will  denote  the  length  of  a  vector 
x,  while  |  Z  |  will  denote  the  determinant  of  a  (square)  matrix  X.  An  asterisk  attached 
to  a  letter  will  indicate  that  the  quantity  in  question  is  a  random  variate. 
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Sets  of  points  will  always  be  Borel  sets,  either  by  hypothesis  or  as  the  result  of  proofs. 
Except  when  the  contrary  is  explicitly  stated,  functions  will  always  be  single  valued.  All 
many-valued  functions  will  be  finitely  many  valued.  All  functions  will  be  Borel  measur¬ 
able  (except  when  explicitly  stated)  either  by  hypothesis  or  as  the  result  of  proofs. 

3.  The  extended  Slutzky-Frechet  theorem 

The  extended  Slutzky-Frechet  theorem  is  the  principal  tool  of  this  paper.  Before  enun¬ 
ciating  it,  we  recall  some  familiar  properties  of  one-valued  random  variates. 

A  random  variate  x*  in  a  space  36  is  defined  by  a  real  nonnegative  countably  additive 
set  function 

(3.1)  F  [?]  =  Pr  { x*  €  ?} 

defined  for  all  Borel  sets  £  of  36  and  satisfying  F[36]  =  1.  F[f]  is  called  the  probability 
set  function  of  x*.  In  the  particular  case  when  j  consists  of  all  points  whose  coordinates 
do  not  exceed  the  corresponding  coordinates  of  a  prescribed  point  x,  we  write  F[j]  = 
F(x),  and  call  F(x)  the  cumulative  distribution  function  of  x*.  Evidently  F[y]  uniquely 
determines  F(x),  and  the  converse  is  a  consequence  of  Lebesgue’s  theory  of  integration. 
Thus  we  can  uniquely  specify  a  random  variate  by  means  of  its  cumulative  distribution 
function.  If  x*,  v  =  1,  2,*  •  *,  is  a  sequence  of  random  variates  in  36  specified  by  the  re¬ 
spective  cumulative  distribution  functions  F,(x),  and  if  lim  Fy(x)  =  F(x)  except  per- 

*— >00 

haps  at  the  hyperplanes  of  discontinuity  of  the  cumulative  distribution  function  F(x), 

we  say  that  x*  converges  in  distribution  to  x*  and  write  dlim  x*  =  x*.  If  y  —  y(x)  is 

*—►00 

a  (one-valued  Borel:measurable)  function  carrying  points  x  of  36  into  points  y  of  another 
space  §),  the  function  of  a  random  variate  y(x*)  is  defined  to  be  the  random  variate  speci¬ 
fied  by  F[i/-1(t))]>  where  y_1(t))  is  the  set  of  all  x  satisfying  y(x)  £  t).  Cramer  [3]  gives 
a  proof  that  this  is  a  consistent  definition,  that  is,  that  F[y_10l)]  is  indeed  a  probability 
set  function  on 

We  now  consider  the  terminology  for  many-valued  functions  of  a  random  variate 
which  will  be  required  in  the  extended  Slutzky-Frechet  theorem. 

An  w-valued  function  y(x)  is  a  set  of  n  points  (not  necessarily  distinct)  in  g)  corre¬ 
sponding  to  each  given  point  x  of  36.  An  indexing  of  y(x)  is  a  system  of  n  one-valued 
functions  yj(x),j  =  1,  2,  •  •  • ,  w,  such  that,  for  every  given  x  of  36,  the  n  points  y(x)  coin¬ 
cide  with  the  n  points  y/x)  with  due  regard  to  multiplicity.  If  there  exists  at  least  one 
such  indexing  in  which  each  yfx)  is  a  Borel-measurable  function,  then  y(x)  is  called  an 
n-valued  Borel-measurable  function ,  and  the  indexing  in  question  is  called  a  Borel-meas¬ 
urable  indexing.  Notice  that  an  w- valued  Borel-measurable  function  may  possess,  besides 
this  Borel-measurable  indexing,  other  indexings  which  are  not  Borel  measurable.  The 
w-valued  function  y(x)  is  said  to  be  continuous  at  the  point  xo  if,  given  any  prescribed 
e  >  0,  we  can  find 

(i)  an  indexing  y,-(x)  of  y(x), 

(ii)  a  number  ij  =  ri  (e,  x©)  >  0,  and 

(iii)  a  permutation  ttx,  7t2, •  •  •,  irn  of  the  integers  1,  2,-  •  •,  w 
such  that 


(3.2) 


ItfyUo)  -y*,.(x)  1  <€ 
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simultaneously  for,;  =  1,  2 whenever  x  satisfies  |x  —  x0|  <  ij.  The  permutation 
in  (iii)  is  allowed  to  depend  on  e,  x0,  and  x.  Then  «-valued  function  y(x)  is  said  to  be 
almost  certainly  continuous  with  respect  to  the  one-valued  random  variate  x*  if  there  exists 
a  Borel  set  jo  satisfying  F[y0]  =  1,  where  F  is  the  probability  set  function  of  x*,  such 
that  y{x)  is  continuous  at  each  x0  £  Jo-  If  yj(x),j  =  1,  2, •••,»,  is  a  Borel-measurable 
indexing  of  an  w-valued  Borel-measurable  function,  and  if  x*  is  a  one-valued  random 
variate  with  a  probability  set  function  F[y],  the  one- valued  random  variate  in  g)  speci- 

n 

fied  by  the  probability  set  function  n~l  is  denoted  by  y(x*)  and  called  the 

j=i 

condensation  of  the  many-valued  random  variate  y(x*).  Evidently  this  definition  is  in¬ 
dependent  of  the  particular  Borel-measurable  indexing  used  therein. 

We  can  now  state  the  extended  Slutzky-Frechet  theorem. 

Theorem  3.1.  If  y(x)  is  a  given  n-valued  Borel-measurable  function  which  is  almost  cer¬ 
tainly  continuous  with  respect  to  a  given  random  variate  x*,  then  dlim  x*  =  x*  implies 

dlim  y(x*)  =  y(x*). 

v-*co 

See  Hammersley  [10]  for  a  proof  of  this  theorem.1 

The  motive  for  condensing  a  many- valued  random  variate  is  to  overlook  its  indexing, 
since  the  indexing  is  an  extraneous  artifice  for  handling  the  different  values.  In  some  re¬ 
spects,  however,  condensation  is  too  drastic  a  way  of  disregarding  the  indexing.  To  ob¬ 
tain  a  less  extreme  concept  consider  any  given  Borel  set  t)  in  g).  For  an  indexing  y/x), 
let  n*  —  n{x*)  denote  the  number  of  values  of  j  such  that  yfx*)  £  t).  Clearly  n*  is 
independent  of  the  particular  indexing  used.  The  moment-generating  set  function  is  defined 
to  be  M(t,  ty)  =  F[exp  (/»*)].  The  greater  flexibility  of  this  concept  is  paid  for  by  in¬ 
creased  mathematical  difficulties;  for,  although  M(t,  ty)  is  a  set  function  of  1),  it  is  not  an 
additive  set  function.  The  mathematical  theory  of  nonadditive  set  functions  is  not  well 
explored. 

The  equation 

(3.3)  M(i,9)  =  [n(  /  +«■/  )]dG(»i,  »„•••,»,) 

'-1  X  Vj  €  V 

provides  a  formal  expression  for  M .  Here  G  is  the  joint  cumulative  distribution  function 

1 1  take  this  opportunity  of  correcting  a  flaw  in  the  proof  of  a  preliminary  lemma  (theorem  1  of  [10]). 
The  flaw  is  concealed  by  the  notation  used  for  the  joint  determination  G,  which  in  fact  depends  on  v,  6, 
and  e.  Instead  of  the  final  relation 

(i)  Gfeo]  >  1  —  e  ,  v  ^  vo  (5,  t)  , 
it  would  have  been  better  to  write 

(ii)  G[3oh  S,  «]  >  1  —  e  ,  v  ^  ><o(5,  e)  . 

What  had  to  be  proved  was  the  existence  of  some  G,  depending  upon  v  but  independent  of  5  and  «,  such 
that 

(iii)  GfeoM  >  1  —  «  ,  v  ^  »-o(5,  «)  . 

The  original  proof  established  (ii)  but  not  (iii).  However  a  reductio  ad  absurdum  argument  suffices  to 
yield  (iii).  For  suppose  the  desired  result  false.  Then  there  exist  fixed  strictly  positive  6,  e  such  that 
(iii)  is  false  for  infinitely  many  v  for  every  joint  determination  G.  Since  5,  e  are  fixed,  (ii)  provides  the  nec¬ 
essary  contradiction. 
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for  the  indexing  yi(x*),  y-z(x*),  •  •  • ,  yn(x*) ;  and  the  product  of  the  integral  signs  is  to  be 
expanded  formally  before  integration.  Thus  for  n  =  2 


(3.4)  M ( t ,  t ))  =  J'jfdGiyi,  y 2)  +  el  f f  dG(yi,  yf) 

vi  e  v 

1/26  t) 

+  el  f  JdG  (yi,  y2)  +  e2t  f  f  dG(yu  yf) . 


€  D' 


y  1  €  9 
1/2  €  t)' 


1/2  €  p 


Indeed,  the  justification  of  (3.3)  is  apparent  as  soon  as  the  formal  expansion  is  made. 
Theorem  3.1  is  a  special  case  of 

Theorem  3.2.  If  y(x)  is  a  given  n-valued  B orel-measurable  function  which  is  almost  cer¬ 
tainly  continuous  with  respect  to  a  random  variate  x*,  if  dlim  x*  =  x*,  and  if  M(t,  ty)  and 

v—>  00 

Mv(t,  X))  are  the  moment-generating  set  functions  of  y(x*)  and  y(x*)  respectively,  then,  for 
every  open  set  X),  lim  Mp(t,  Xf)  =  M(t,  X)). 

v—*oo 

This  theorem  follows  without  difficulty  from  the  result  established  by  Hammersley 
(see  1.  10,  p.  256  in  [10]).  The  requirement  that  i)  should  be  an  open  set  can  be  relaxed 
somewhat;  roughly  speaking,  all  that  is  needed  is  that  M(t,  fy)  should  not  be  a  discon¬ 
tinuous  function  of  l)  at  ty,  but  a  precise  formulation  of  this  idea  is  rather  cumbrous. 

Evidently  the  probability  set  function  of  y(x*)  is  the  coefficient  of  t  in  the  expansion 
of  M(t,  t i)/n  in  powers  of  t. 


4.  Borel  measurability  and  continuity  of  the  zeros  of  a  polynomial 

In  this  section  we  prove  that  the  zeros  of  a  polynomial  are  Borel  measurable  and  (with 
a  single  exception)  continuous  functions  of  the  coefficients  of  the  polynomial. 

Lemma  4.1.  There  exists  a  fixed  indexing  Zi,  z2,  •  •  * ,  zm  of  any  given  fixed  set  of  m  distinct 
points  in  the  complex  plane  such  that 

(4.1)  |  Zi  -  eri /*  |  >  |  z2  —  eTi/*  \  >■■■>  \  zm-  eri^  \  >  0 

for  all  sufficiently  large  positive  integers  s,  say  s  ^  So(zi,  z2,‘  •  • ,  zm). 

In  enunciating  this  lemma  we  make  no  attempt  to  claim  that  S0  is  a  Borel-measurable 
function  of  Z\,  z2,  •  •  • ,  zm.  The  reason  for  this  remark  will  appear  during  the  proof  of  theo¬ 
rem  4.1  below. 

Lemma  4.1  effectively  states  that  an  ordering  can  be  set  up  for  the  m  points,  and  it  is 
therefore  enough  to  consider  the  case  m  =  2.  The  result  is  obvious  when  Z\  and  z2  are  at 
different  distances  from  z  =  1.  When  j  zx  —  1 1  =  |  z2  —  1 1 ,  zi  is  the  name  given  to  the 
point  with  the  smaller  imaginary  part,  except  that,  when  the  two  imaginary  parts  are 
equal,  Z\  is  the  name  of  the  point  with  the  larger  real  part. 

Theorem  4.1.  There  exists  an  indexing  Z\,  z2,  •  •  • ,  zn  of  the  zeros  of  the  polynomial 

n 

^  CjZj,  where  c3-  =  a}-  +  ibj,  such  that, for  each  k  =  1,  2, •••,»,  z*  =  {xk,  yk]  is  a  one- 
r=  0 

valued  Borel-measurable  function  of  c  =  {a0,  ai,  -  •  -  ,  an,  bo,  bi,-  •  •,  bn).  Moreover,  for  this 
indexing  |  zk  —  1 1  is  a  nonincreasing  function  of  k. 

First  and  until  further  notice  assume  that  c  is  fixed  and  c„  =  1.  Then  the  zeros  of 

n 

^  CjZ]'  are  fixed  and  all  finite.  Suppose  that  there  are  m  distinct  zeros  amongst  them, 
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and  use  the  indexing  of  lemma  4.1  to  name  these  distinct  zeros  Zi,  z2,  •  •  • ,  zm.  Since  zi,  z2, 
•  •  • ,  zm  are  uniquely  determined  by  c,  the  function  50  of  lemma  4.1  is  a  function  of  c,  say 
50(zi,  z2,  •  •  • ,  zm)  =  5(c).  We  have  not  yet  proved  that  zh  z2,  •  •  • ,  zm  are  Borel-measurable 
functions  of  c,  so  that  we  cannot  and  do  not  assert  that  5(c)  is  a  Borel-measurable  func¬ 
tion  of  c.  This  explains  why  we  did  not  bother  to  prove  that  So  was  a  Borel-measurable 
function  of  Zi,  z2,  •  •  • ,  zm. 

For  each  integer  s  =  1,  2,  -  •  •  define  numbers  C]S  by  means  of  the  identity  in  z 
(4.2)  V  £,(«  +  «•!/•)>=  ]£Cy,z>'. 

7  =  0  7=0 


For  t  =  1,  2,  •  •  • ,  n  let  ut  be  an  arbitrary  real  number  satisfying  0  ^  ut  ^  1.  For  s,  t  = 
1,  2,  •  •  •  define  Ztt  as  a  function  of  u  —  («i,  «2,  •  •  • ,  un\  by  means  of  the  recurrence  rela¬ 
tions 

f  ut  for  n , 


(4.3) 


Then  define 
(4.4) 


Z't  =  { 


n— 1 

-2 

7-0 


Cj'Z 


t,  t—n+j 


for  /  >  n  . 


tf.tH  z“ 


for  Ztt  0  , 


Finally  write 
(4.5)  /(c) 


0  for  Zat  =  0  . 


=  lim  f  f  ■  f  lim  ( eri^  +  Utt)  duxdu2  •  •  •  dun . 

a — >co  J o  Jo  Jo  <— »oo 


We  shall  prove  that /(c)  exists,  is  a  Borel-measurable  function  of  c  subject  to  cn  —  1, 
and  equals  Z\. 

By  (4.2),  cn  —  1  implies  Cnt  =  1.  Hence  the  linear  difference  equation  (4.3)  has  the 
solution 

m 

(4.6)  Z.^^PMiZk-  <:"'•)• ,  t=  1,2,  -, 

k  =  l 


where  Pk(t)  is  a  polynomial  in  t  whose  coefficients  are  (not  necessarily  Borel  measurable) 
functions  of  u,  for  the  distinct  roots  of  the  right-hand  side  of  (4.2)  are  z  =  zk  —  exp 
( ici/s ),  k  =  1,  2,-  •  •,  m.  The  polynomials  Pk(t)  have  degree  less  than  n  and  are  all 
identically  zero  if  and  only  if  u  =  0.  When  u  7^  0,  let  K  be  the  smallest  value  of  k  such 
that  Pk(t )  is  not  identically  zero.  Consider  any  fixed  value  of  s  ^  5(c).  By  lemma  4.1 

(4.7)  |  zk  -  e™/*  |  >  |  z*+1  -  |  >  •  •  •  >  |  zm  -  |  >  0 


Hence,  from  (4.6), 

(zK—  eTi/*  for  u  0  , 

(4.8)  lim  =  < 

t_>°o  (0  for  u  =  0  . 


Now  K  is  a  function  of  u,  but  we  show  that  K  —  1  for  almost  all  u.  From  (4.3)  and  (4.6) 


ut 


m 


^Pi(0(Zi-e',/')‘, 

fc-1 


(4.9) 


t  =  1,  2,  •  ■  •  ,  n  . 
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There  are  altogether  n  coefficients  in  all  the  polynomials  P*  taken  together.  If  any  one 
of  these  coefficients  is  zero,  (4.9)  shows  that  there  is  a  linear  relation  between  ui,  u*,  ■  •  • , 
un;  that  is,  u  lies  on  a  hyperplane.  Hence,  except  perhaps  when  u  lies  on  one  or  more  of 
n  hyperplanes  in  the  unit  hypercube,  none  of  the  polynomials  P*  is  zero.  Hence, 
from  (4.8), 

(4.10)  nr  lim  (  U3t -h  e1"7*)  du\du%  •  •  •  dun  =  Zi ,  s  ^  S(c) . 


By  (4.5), /(c)  exists  and/(c)  =  Zi  for  each  fixed  c  subject  to  cn  =  1.  It  remains  to  prove 
that  /(c)  is  a  Borel-measurable  function  of  c.  We  therefore  allow  c  to  vary  subject  to 
the  restriction  c„  =  1.  Evidently,  by  (4.2),  the  Cja  are  Borel-measurable  functions  of  c 
and  s;  hence,  by  (4.3)  and  (4.4),  U,t  is  a  Borel-measurable  function  of  c,  s,  and  u.  The 
result  follows  from  (4.5),  since  the  limit  and  the  integral  of  Borel-measurable  functions 
are  Borel  measurable. 

Next  we  remove  the  restriction  cn  =  1.  If  cn  0,  we  may  divide  all  the  coefficients  of 

W 

CjZi  by  cn  without  affecting  z\.  If  cn  =  0,  we  define  Zi  =  { + 00 ,  +  °°  } .  It  follows  that 

J-0 

zi  =  yi\  is  a  Borel-measurable  function  of  c  without  restriction  upon  c. 

n— 1 

For  the  next  zero,  we  repeat  the  foregoing  process  for  the  polynomial  7 ,zj,  whose 

J-0 

n 

zeros  are  those  of  c]  zi  excepting  one  zero  at  z  =  z\.  The  yj  are  defined  by 

J-0 


(4.11) 


J  0  for  j  —  n  , 

l  Cj+i  -  Zi7j+i  for  /  =  n  -  1,  n  -  2,  •••  ,  1,  0  , 


provided  c„  5^  0,  and  by  7/  =  Cj,  j  =  0,  1,  •  ■  • ,  n  —  1,  if  cn  =  0.  Since  a  Borel-measur¬ 
able  function  of  a  Borel-measurable  function  is  Borel  measurable,  the  7 ,•  are  Borel- 

W 

measurable  functions  of  c.  Hence  the  next  zero  of  cjz1'  is  a  Borel-measurable  func- 

J-0 

tion  of  c. 

Further  repetition  of  this  process  provides  the  required  indexing  and  shows  that  the  z* 
are  Borel-measurable  functions  of  c;  the  method  of  procedure  allied  with  lemma  4.1 
shows  that  \zk  —  1  ]  is  a  nonincreasing  function  of  k. 

Lemma  4.2.  Let  nbe  a  given  integer  greater  than  zero,  and  let  R  and  ci  be  given  real  num¬ 
bers  satisfying  0<d^2  ^  R  <  «.  Define 


(4.12) 

Let 


i  =  1,  2,  ••  •  ,  n. 


(4.13)  Pn(z)  =  ^PnjZ’  =  0  ,  Qn(z)  =  ^2  qnjz’=  0,  pnn  =  qnn=  l 

j'-o 


J-0 


be  two  algebraic  equations  {with  complex  coefficients )  having  no  roots  outside  the  circle  \z\  = 
R.  Then,  if 


(4.14) 


I  Pnj  —  ?nj  |  ^  «n+l  > 


/  =  0,  1,  ••  •,«  -  1 
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it  is  possible  to  arrange  the  roots  a  of  Pn(z )  =  0  and  the  roots  P  of  Qn(z)  =  0  into  some  set 
of  n  pairs  (ay,  Pj)  such  that 

(4.1 5)  |  ay  ft,  |  <  «y  >  j =  1  >  2,  •  •  •  ,  n  . 

Since  \  ^  «i/2 R  =  >  0  it  is  clear  from  (4.12)  that  ^  ^  €2  >  e3  >  •  •  ■  >  en+i  >  0; 

hence 

(4.16)  2ey+1^2^1(2ei+1)1/J',  i=l,  2,  •••,». 

If  X  is  any  root  of  Pn(z )  =  0  or  of  Qn(z)  =  0  we  have  |  X  |  ^  R,  and  hence 

(4.17)  S  <W 

kmmQ  Jc—Q  A  ^ 

since  JR  ^  2. 

Let  j8n  be  any  specified  root  of  Qn(z)  =  0.  Then 

(4.i8i  n  («*-•»)  -  !*•<«  1  -  -<?.<«  1  - 1  s  ^ 

j— 1  j— 0 

|fl[|<2*MR" 

0 

by  virtue  of  (4.14)  and  (4.17).  It  follows  that  there  is  at  least  one  root  of  Pn(z)  =  0, 
which  we  denote  by  an,  such  that 

(4.19)  |/Sn-an|  <  (2 en+1)^R. 

Now  define  the  polynomials  Pn-i(z)  and  Qn-i(z)  by  the  identities 

n— 1 

pn  (  z)  SS  (  Z  ~  an)Pn—i  (  z)  =  (  z  -  an)  2  P»-l.  i  Z’  >  £»-l,  n-X  =  1  , 

y— 0 

(4.20) 

n— 1 

Q»(  z)  =  (  Z  -  /3n)  @B-i  (  z)  ss  (  Z  -  /S„)  2  0»-l.  i  »  ?n-l,  n-1  =  1  • 

7“0 

The  roots  of  JPn-i(z )  =  0  and  Qn-i(z)  =  0,  being  a  subset  of  the  roots  of  P„(z)  =  0  and 
Qn(z)  =  0,  do  not  lie  outside  the  circle  |  z  j  =  R.  Upon  identification  of  the  coefficients 
in  (4.20), 

(4.21)  (pn—it  j  1,  y)  =  ( pn  y+i  Qn,  y+i)  "4*  &n(pn— 1>  y+i  9»— 1  y+i) 

"4"  (®»”  Pn)  ?n-l,  y+1  • 

Now  ±  5n_i,  y+i  is  the  sum  of  the  products  of  the  roots  of  Qn-i(z)  —  0  taken  n  —  j  —  2 
at  a  time,  so 

(4.22)  |  | 

Take  the  modulus  of  (4.21)  and  employ  (4.14),  (4.19),  and  (4.22)  to  yield 

(4.23)  |  pn—i,  y  ?»—  1,  y  |  ^®»+i~HP  I  Pn—  1  y+i  1,  y+il 

+  (  2«n+l)  l/nR  (y  ^  j)  . 
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In  this  inequality  write  j  +  k  for  j,  multiply  through  by  Rk,  and  sum  over  0  ^  k  ^ 
n  —  j  —  2.  This  yields 

(4.24)  ^.R*  |  Qn  1  f  j+Jc  |  ^  |  Pn — lt  j+fc+l  Qn — 1,  j+k+l  | 

A:— 0  A;=0 

+"§  k**+o+*+  ,)(2 . 

the  first  sum  on  the  right-hand  side  having  the  upper  limit  n  —  j  —  3  instead  of  n  — 
j  —  2  since  pn-i,n-i  =  qn-i,n-i  =  1-  Cancelling  like  terms  from  each  side  of  the  in¬ 
equality,  we  find 

(4.25)  <  V’  j «»+i 2?*  +  ( 26,.+!) '/" J 

<£n+l  £**  +  (2 «„+,)  ■/.*»-.  W*  ,  *) 

A:=0  \  &  / 

<  2e„+1R“-’  +  (2e„+1)‘/"(22?)"-i 

5  2(26n+1)l/”(2J?)”-1g  (2JJ)%^1  =  «., 

upon  using  (4.17),  (4.16),  and  (4.12).  But  (4.25)  is  simply  (4.14)  with  n  —  1  for  n.  Hence 
we  may  repeat  the  process  successively  until,  for  any  given  order  of  the  roots  /3y,  we  have 
picked  out  the  roots  ay  in  such  an  order  that 

(4.26)  I  ftj  ~  ay  |  <  (2cy+i)1/^<€y 

by  virtue  of  (4.19),  (4.12),  and  R  ^  2. 

The  proof  given  above  follows  the  general  lines  of  an  earlier  attack  by  Coolidge  [2], 
but  care  has  been  taken  to  ensure  the  pairing  of  the  roots.  For  our  purposes  it  would  not 
have  sufficed  to  show  that  every  root  of  P„(z)  =  0  is  near  some  root  of  Qn(z)  =  0  and 


vice  versa. 


Lemma  4.3.  None  of  the  roots  of  the  equation  g(z )  =  ^  cjzj  =  0  lie  outside  the  circle 

y=o 

|  z  |  =  1  4-  M/L,  where  L  —  \  cn  |  and  M  =  max  |  cy  | . 

0^7<n 

For  suppose,  on  the  contrary,  that  g(zo )  =  0  and  H  =  |zo|  >  1  4*  M/L  ^  1.  Then 

I  n  n— 1 

2  CA  ^LHn-M'£Hi  =  LHn-M(Hn-l)  /  (H -1)  .SinceH>  1,  we 

7=0  7=0 

deduce  0  2s  L{H  —  1)  —  Mil  —  H~n)  ^  L(H  —  1)  —  M,  and  thence  the  contradic¬ 
tion  H  ^  1  +  M/L. 

n 

Theorem  4.2.  The  n-valued  function  consisting  of  the  zeros  of  the  polynomial  c]  zi  —  0 

7  =  0 

is  continuous  at  any  point  c  such  that  cn  9^  0  and  cn- 1,  cn- 2,  •  •  • ,  Co  are  finite. 

Suppose  e  >  0  is  prescribed.  Consider  any  fixed  c  such  that  cn  5*  0  and  cn- 1,  c„_ 2,  •  •  • , 
Co  are  finite.  Choose  5  >  0  so  that 


(4.27) 


e  >  5  >  0  ,  Slc>i=T’ 
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Define 

(4.28) 


(4.29) 


*(*)"“  2^**1  h(z)  =  '£djz’,  \c-d\<v . 


(4.30) 


cj  —  dj  |  <  v  • 


In  view  of  the  definition  of  continuity  in  section  3,  it  is  enough  to  show  that  the  zeros  of 
g(z)  and  h(z)  can  be  arranged  in  pairs  such  that  the  distance  in  the  complex  plane  be¬ 
tween  the  two  members  of  any  pair  is  less  than  8. 

Lemma  4.3,  in  conjunction  with  (4.27),  (4.28),  and  (4.29),  shows  that  none  of  the 
roots  of  g(z )  =  0  and  h(z)  =  0  lie  outside  the  circle  |z|  =  1  +  (5_1  +  17)/ (5  —  17)  < 
2 /S2.  In  the  notation  of  lemma  4.2,  we  take  R  =  2/8 2  and  €1  =  8.  This  ensures  0  <  €1  ^ 
2  5s  R  <  00 .  The  roots  of  g(z)  =  0  and  h(z)  =  0  lie  inside  the  circle  \z\  =  R.  Write 
pnj  =  Cj/cn  and  qnj  =  dj/dn;  so  that  pnn  =  qnn  =  1,  and  Pn{z)  =  0  has  the  same  roots 
as  g(z)  =  0  while  Qn(z)  =  0  has  the  same  roots  as  h(z )  =  0.  Now,  for  0  ^  j  ^  n  —  1, 

/.n..  1  .  1  Cj[cn~ 4"  ( dn  cn )  ]  cn[  Cj  ( dj  c/)  1 

(4.31)  - c-U»--rU;--cn)) - 


Pnj  qnj  |  — 


Cj  ( dn  Cn)  . 

=  Cn[cn+(dn- cn)  ]  ^ 

^  S^y  ■  V  g2q 
~~  8  (  8  —  ij)  8  —  y~  83 


dj  ci 
cn+  ( dn  -  Cn ) 


by  virtue  of  (4.30)  and  (4.27)  and  8  ^  Finally,  from  (4.12), 

(4.32)  log  €j+i  =  j  log  tj  -  j 2  log  ( 2 R) . 

Multiply  this  equation  by  n\/j\  and  sum  over  1  fg  j  n,  to  yield 

(4.33)  log  [log «,-{ log  ( 21?) } 

00  • 

^  »!  [log  Cl  -  {log  (27?)  }  2  ;]  =  nl  Hog  Ci  —  2e  log  (2 R)  ] 

6  »!  log  =  log  (^?)  • 

The  result  now  follows  from  (4.31),  (4.33),  and  lemma  4.2,  since  the  tj  form  a  decreas¬ 
ing  sequence. 

5.  Real  isomorphs  of  complex  matrices 

Any  complex  matrix  Z  can  be  expressed  in  the  form  Z  =  X  +  iY,  where  X  and  Y 
are  real  matrices.  We  define  the  real  isomorph  of  Z  to  be 

'-<?  -D- 


(5.1) 


IOO 
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Wedderburn  (see  p.  100  in  [18])  uses  this  device  in  the  study  of  Hermitian  matrices:  a 
matrix  is  Hermitian  if  and  only  if  its  real  isomorph  is  symmetric.  We  use  the  word  “iso- 
morph”  in  view  of  the  easily  verified  relations 


(5.2) 


Zx-\-Z2  =  Zx-\-Z2  , 
Z \Z 2  —  Z \Z 2  , 


which  are  valid  in  the  sense  that,  when  one  side  of  one  of  the  equations  (5.2)  exists,  the 
other  side  of  that  equation  exists  and  the  two  sides  are  equal.  The  matrix  identity 


(5.3) 


(I  iI\(X  ~Y\(I  -iI\  =  /X+iY  0  \ 
VO  /Ar  lAo  i)  \  Y  X-iYj 


implies 

Theorem  5.1.  If  Z  is  a  square  matrix ,  the  determinant  of  2  equals  the  square  of  the 
modulus  of  the  determinant  of  Z. 

We  have  already  mentioned  that  the  complex  number  z  —  x  +  iy  may  be  regarded 
as  a  vector  z  =  {x,  y} .  Corresponding  to  the  latter  form,  we  shall  write  dz  as  an  abbre¬ 
viation  for  dxdy,  while  dz  will  represent  the  ordinary  complex  differential  element.  If  a 
multiple  integral  in  a  number  of  complex  variables  w i,  w2,  •  •  • ,  wm  is  to  be  transformed 
to  an  integral  in  the  complex  variables  Z\,  z2,  •  •  • ,  zm  by  means  of  the  relations 


(5.4) 


Wj=Wj(z1,  z 2,  •••  ,  Zm), 


j  1>2 ,  *  ■  ■  >  m  i 


the  appropriate  Jacobian  d(wx, •  ■  •,  wm)/d(zx,-  •  •,  zm)  will  consist,  in  the  usual  way,  of 
the  determinant  whose  elements  are  dWj/dzk.  Let  us  suppose  that,  for  each  J  and  k,  w,  is 
an  analytic  function  of  zk  when  Zi,  •  •  • ,  zk~ i,  zk+x,  •  •  • ,  zm  are  held  fixed.  The  familiar 
Cauchy-Riemann  equations  assert,  for  an  analytic  function  w(z)  =  u(z)  -f-  iv(z),  that 

__  du_dv_^dw  d  v  _  _  du  _  dw 

{'}  Jx~dy~Ji  Jz'  Jx~  dy~  ^Tz' 


where  and  3  denote  the  real  and  imaginary  parts  of  the  ensuing  expression.  From  (5.5) 
and  (5.1),  it  follows  that,  if  J  is  the  Jacobian  matrix  of  the  transformation  from  dwx  •  •  • 
dwm  to  dz\  ■  •  -dzm,  then  J  is  the  Jacobian  matrix  of  the  corresponding  transformation 
from  dw\  ■  •  • dwm  to  dzv  *  • dzm •  Hence,  theorem  5.1  implies 
Theorem  5.2.  If 

(5,6)  dwx-'-  dwm  =  d.  dzi  •••  dzm, 

o  t  z i,  •  •  •  ,  zm) 


then 

(5.7) 


dwx  •  •  •  dwm  = 


d  ( wx ,  •  • 

*  )  Wm) 

d  (Zi,  •• 

)  zm) 

dzi  *  dzm . 


Further  than  this,  the  isomorphism  exhibited  in  (5.2)  enables  us  to  construct  directly 
the  analogues  of  all  the  familiar  theorems  on  Jacobians,  which  rest  on  the  ordinary  prop¬ 
erties  of  matrix  multiplication.  For  instance, 

Theorem  5.3.  If  the  variables  wx,  ■  •  • ,  wm  are  related  to  the  variables  zx,  •  *  • ,  zm  by  analytic 
equations 

(5.8)  Fj(wx,  •  •  •  ,w„,  zx,  •  •  •  ,  zm)  =  0  ,  J=l,  2,  •••,*», 


then 

(5.9)  dwx 


d  (Fx,  •  •  •  ,  Fm)  /d{Fx,  •••  ,Fm)  * 
d(zx,  ,  zm)  /  d(wx,  •••  ,wm ) 


dzx  "  d z m 
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6.  The  moment-generating  set  function  of  the  zeros  of  a  polynomial  with  general¬ 
ly  distributed  coefficients 

We  derive  a  formula  for  the  moment-generating  set  function  of  the  zeros  of  a  poly- 

n 

nomial  g*(z)  =  ^  c*z]\  when  the  joint  cumulative  distribution  function  of  the  coeffi- 
;'-0 

cients  has  the  arbitrary  form 

(6.1)  W{c)  =W(a0,  «!,•••,  fl„,  b0,  bi,  ••• ,  bn) . 


We  first  obtain  this  moment-generating  set  function  when  certain  restrictions  are  placed 
on  W,  and  then  later  invoke  theorems  3.1  and  3.2  to  relax  these  restrictions.  Accordingly 
we  suppose  until  further  notice  that  W  possesses  a  continuous  frequency  function 


(6.2) 


w(c ) 


_ din+*W(c) _ 

daedal  •  *  •  ddnd  bodbi  •  •  •  d  bn  ’ 


and  that  all  moments  of  c*  exist. 

Let  Gt  denote  the  (2 n  +  2)-dimensional  Euclidean  space  consisting  of  points  c,  and  3 
the  complex  plane.  For  p  =  1,  2,  •  •  • ,  r  let  dzp  be  the  small  rectangle 

(6.3)  xp<  x  ^  xp+dxp  ,  yP<y^yP+dyp, 

where  zp  =  xp  +  iyp  are  prescribed.  We  are  going  to  evaluate  the  probability  that 
g*(z)  =  0  has  precisely  one  root  in  each  of  the  small  rectangles  dzp.  The  probability  will 
be  written  in  the  form 

(6.4)  P  (Zi»  •  *  •  >  Zr)  dll  ’  ’  •  dzr  . 

This  expression  corresponds  to  Ramakrishnan’s  product  density  [16],  introduced  in  the 
treatment  of  cosmic  ray  phenomena.  However,  in  Ramakrishnan’s  work  there  was  a 
certain  degree  of  independence  between  behaviours  of  distinct  quanta;  whereas  in  the 
present  context  we  do  not  enjoy  any  such  simplification  because,  when  one  root  of  g*  is 
known,  the  conditional  distribution  of  the  remaining  roots  is  affected. 

For  fixed  z  =  x-\-  iy,  the  equations 

(6.5)  ftg(z)  =  3g(z)  =0 

define  a  2«-flat  f*  in  (S.  As  z  varies,  the  flats  f,  develop  a  twisted  regulus  r.  The  generator 
of  r  lying  in  fz  is  a  (2 n  —  2)-flat  g*  with  equation 

(6.6)  $g(z)  =  Sg(z)  =  Wg'(z)  =  Sg'(z)  =  0  . 

The  Cauchy-Riemann  equations  provide  the  reason  why  g*  is  a  (2 n  —  2)-flat  and  not 
a  (2 n  —  4)-flat,  as  one  might  at  first  sight  imagine. 

If  g*(z)  =  0  has  at  least  one  root  in  each  of  the  rectangles  dzp,  p  =  1,  2,  •  •  • ,  r,  then  c* 
must  lie  at  the  intersection  of  the  r  2«-flats  f*p.  If,  further,  one  of  the  rectangles  dzp 
(say  dzx)  contains  more  than  one  root  of  g*(z)  =  0,  then  c*  must  also  lie  on  g^.  Since 
w(c)  is  supposed  continuous,  it  follows  that  the  probability  that  each  of  the  dzp  con¬ 
tains  precisely  one  root  differs  by  terms  of  higher  order  than  dzx  •  •  •  dz,.  Hence 

(6.7)  P(zi,  •  •  •  ,  zr)  dzi  •  ■  •  dzr  =  J^w(c)  dd0  ■  •  •  ddndb0  •  •  •  dbn, 


where  c  is  the  intersection  of  the  f,,  for  zp  £  dzp.  To  evaluate  (6.7),  we  transform  from 
the  coordinate  system  (ao,  di,  •  •  • ,  dn,  bo,  b\,  •  •  • ,  bn)  to  the  coordinate  system  (xi,  xi,  •  •  • , 
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xr,  aT,  ar+i,  •  •  -  ,  an,  yi,  >v  •  •,  yr,  br,  bT+i,  •  •  •,  bn).  The  relationships  between  these  sys¬ 
tems  are 

(6.8)  g(zp )  =0  ,  p=  1,  2,  •  ••  ,  r  . 

Now  this  is  the  situation  envisaged  by  theorem  5.3.  The  appropriate  Jacobian  is  accord¬ 
ingly 


3(*.  •• 

g)  / 

a (g,  •  - 

• .  g) 

2 

d  (zi,  •  • 

•  ,  *r)  / 

d  ( Co.  *  *  * 

g' (2,) 

0 

>  cr— i) 

0 

*'U2) 

...  o 

...  0 

/ 

1  z  •  •  •  zr— 1 

i  i 

1  v  • 

0 

0 

•••  f'Ur) 

/ 

1  «r 

=  n  u'(2P)  i2/  n  i  z*!2- 

P= 1  '  1 gp<?gr 

In  (6.9)  the  denominator  is  to  be  taken  as  unity  if  r  =  1.  Thus 


(6.10)  P(i„  •■•,*,)  -  / 

o(zp)=0;  p  =  l, 


n  i  s' Up)  i2 

— - - w(c)  dcT  •  •  ■  dcn . 

r  J~J  I  ZP  —  zq  I  2 
lgp< q£r 


In  (6.10)  the  conditions  defining  the  range  of  integration  are  to  be  used  to  express  c0,  ’ '  ’ , 
cr-i  in  terms  of  cr,  •  •  • ,  cn  and  zh  ■  •  ■ ,  zn.  The  integral  (6.10)  exists  because  all  the  moments 
of  w(c )  are  supposed  to  exist. 

We  now  consider  the  important  special  case  r  =  n.  The  integration  in  (6.10)  will  be 
with  respect  to  dcn  only.  Since 


(6.11)  g(z)  =  c„f"J  (z  -  Zj) , 

j=i 

we  have 

(6.12)  «'(«,)  =  (%-*,); 

jVp 

and  hence 


(6.13)  n  i i*i-  n  i  **— *«i4- 

P=1  l*£P<Sg»> 

Thus 

(6.14)  P  ( Zi,  •  •  •  j  zn)  =  T  ]~J  I  Zp“  z«l2l/l  c»|2nw(c) 

Llgp<«^n  J 


In  (6.14)  it  is  supposed  that  the  argument  c  of  w(c)  is  expressed  in  the  form 
(6.15)  an-r  =  9?  (  cnHr)  ,  bn-r  =  3  (  cnHr) ,  r  =  1,2, 
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where  Hr  is  the  sum  of  the  products  of  —  zi,  —  22,  •  •  •,  —  z„  taken  r  at  a  time.  Since  there 
are  n\  permutations  of  Zi,  Z2,  •  •  • ,  zn,  we  have  on  combination  of  (3.3)  and  (6.14) 


(6.16)  M(U)=l[nj  /  +e‘fl]p(.Zl,  ■■■,z.)dzl- 

}=1  Zj  6  V  Zj  e  i 

Another  important  special  case  of  (6.10)  is  r  =  1.  We  obtain 

(6.17)  P(zi)=  J  |  g'(zi)  \2w(c)dci  •••  dcn 

0(zi)  =0 


dZn. 


=  f  \g'\2^(g>  g')dg', 

g(z  1)  =0 

where  \f/(g,  g')  is  the  joint  frequency  function  of  the  real  and  imaginary  parts  of  g*(zi) 
and  g*'(zi).  However  P(zi)dzi  is  the  expected  number  of  roots  in  the  small  rectangle  dZ\\ 
and  therefore  the  condensed  probability  set  function  of  the  zeros  of  g*(z)  is 

(6.18)  F(3)  =—  Cp (zi)  izi. 

n 

Finally  we  can  remove  the  restrictions  originally  placed  on  c*.  For  we  can  find  a  se¬ 
quence  of  random  variables  c*,v=  1,  2,  •  •  • ,  such  that  each  c*  has  moments  of  all  order 
and  a  continuous  frequency  function,  and  such  that  dlim  c*  =  c*.  Corresponding  to 

CO 

ct  we  can  find  M„(t,  3)  or  F,(i)  as  in  (6.16)  or  (6.18);  and,  when  3  isopen, 

(6.19)  =  lim  M  (1,  3) »  F(3)  =  lim  F  (8) 

V — >00  v — ^oo  y 

which  provides  the  required  set  functions  corresponding  to  c*,  in  accordance  with  theo¬ 
rems  3.1  and  3.2.  In  removing  the  restrictions  on  c*,  we  ought  to  assume  (in  order  to  avail 
ourselves  of  theorems  3.2  and  4.2)  that  c*  =  0  has  not  a  positive  probability.  But  if 
c*  =  0  has  a  positive  probability,  there  is  the  same  positive  probability  of  an  infinite 
root  of  g*  =  0;  and  we  can  work  with  the  conditional  probability  given  that  c*  ^  0. 


7.  The  characteristic  functional 

The  characteristic  functional 

(7.1)  C[0]  =£["exp 

offers  a  more  general  approach  than  the  moment-generating  set  function  and  the  method 
of  product  densities  considered  in  section  6.  On  the  other  hand  it  seems  much  less  trac¬ 
table,  and  I  have  been  unable  to  derive  any  useful  results  from  it;  but  I  present  it  never¬ 
theless  in  the  hope  that  the  reader  may  manage  to  manipulate  it  successfully.  In  (7.1), 
6(z )  is  an  arbitrary  function  of  z,  and  z*  are  the  zeros  of  g*(z)  =  0.  We  have  to  express 

n 

^  6(z*)  in  terms  of  g*,  and  then  take  expected  values  in  (7.1)  with  respect  to  the  dis- 

7-1 

tribution  of  g*.  The  first  part  of  this  programme  can  be  achieved,  at  least  symbolically, 
when  6  is  an  integral  function  by  using  Cauchy’s  theorem: 

(^•2)  C  [01  =E  [  lim  exp  {  {2tt)  -1  $  6  (  z)  d  log  g*  ( z)  }  ]  . 

R— |z|=i? 
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The  second  part  of  the  programme  is  however  harder,  since  g*(z)  regarded  as  a  function 
of  z  is  a  deterministic  random  function  whose  value  is  known  for  all  z  as  soon  as  it  is 
known  for  ti  distinct  values  ft,  ■  •  • ,  ft.  Thus,  by  Lagrange’s  formula, 


(7.3) 


g*(z) 


=£«*(?>)  n 

l  A*;' 


2  ~  ft 
ft  -  ft 


can  be  inserted  in  (7.2)  and  the  expected  value  taken  over  the  n-dimensional  joint  dis¬ 
tribution  of  g*(ft)  for  convenient  fixed  ft.  But  the  manipulative  difficulties  seem  severe. 

The  characteristic  functional  reduces  to  the  moment-generating  set  function  for  a  set 
3  if  9  is  taken  equal  to  —  it  times  the  indicator  function  of  3.  If  3  is  bounded  by  a  Jordan 
curve  and  A3  log  g*(z )  denotes  the  variation  of  log  g*(z)  in  describing  this  curve,  we 
may  write 

(7.4)  M (t,  3)  =  -E^exp  Aj  log  g*  (z)  |]  , 

into  which  (7.3)  may  be  substituted  as  before. 


8.  Normally  distributed  complex  coefficients 

We  shall  determine  the  condensed  distribution  of  the  zeros  of  g*(z)  when  the  coeffi¬ 
cients  c*  are  normally  distributed  about  a  mean  y  with  nondegenerate  variance-covari¬ 
ance  matrix  V.  The  frequency  function  of  c*  is  thus 

(8.1)  ( 2  it)  n  *  |  V  ] — exp 


We  require 

Lemma  8.1.  Assuming  conformable  partitioning  and  the  existence  of  the  relevant  inverses, 
and  provided  A  is  symmetric, 


(8.2) 

iaV1  (Oi\ 

A22J 

=  aldii1ai+  (<12  - 

”^21-4ll1fll)  *  (^22  -421-dll  An)  (<*2 

Let 

(8.3) 

II 

t!)  (b1  bO=/’ 

A22/  \0  21  -022/ 

so  that 

(8.4) 

/  AnBn  ~\~AnB2i 
\42i5ii  +^22^21 

illl5l2+^12-B22>|_//  0\ 

A2\B\2~\~ A22B22)  VO  1/  . 

The  upper  right-hand  comer  of  (8.4)  gives 

(8.5)  Bn  —  ~AiiAuBz2, 


and  substitution  of  (8.5)  into  the  lower  right-hand  corner  of  (8.4)  gives 


(8.6)  B22  —  (^22  A21A11A12) 

Similarly  the  upper  right-hand  comer  of  (8.4)  gives  An  =  —  AuBnB^f,  which 
substituted  into  the  upper  left-hand  corner  of  (8.4)  yields 

(8.7)  .An1==2?n  B12B22B21  • 

Since  A  is  symmetric,  so  is  B;  and  (8.5)  gives 

(8.8)  ^uAu.1  =  (^u^tt)  ^ =  (  B12B22  )  =  B22B21  . 
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From  (8.2),  (8.6),  (8.7),  (8.8) 

(8.9)  +  (fl2  — i421^111Al)  f  (  A22  —  AnAuAu)  1  ( fl2  A2lAn1ffi) 

=  Cl  (2?n  —  Bit  Ba1  B21)  fll+  (fl2+522152iai)  ^22  (fl2-l-'ff22152lOl) 

=  a!j?nai  +  2  ai5i2a2  +  02^2202  =  a*Ba  =  <tA~xa  . 

Lemma  8.2.  |A|  —  |An|  |A22  —  A2iAui*Ai2|  • 

This  follows  immediately  from  the  identity 


(8.10) 


—  ^21^111 


A22  ~  A2iAh .4 12 


Let  z  —  x  +  iy  be  a  fixed  complex  number,  and  define 

(Zo-  (1,  ,  zn) 


(8.11) 


Zi  =  (0,  1,  2  z,  •  •  •  ,  «2n_1) 


It  is  easy  to  verify  that  the  column  vector  q  defined  by 

(8.12)  <7=  {&g,  SR*',  3*,  Sg'}  =  Zc  ; 

and  hence  the  real  and  imaginary  parts  of  g*(z)  and  g*'(z)  are  normally  distributed  about 
a  mean  2y  with  variance-covariance  matrix  ZVft.  Here  and  elsewhere  ft  denotes  the 
transpose  of  2  and  not  the  real  isomorph  of  /t.  From  (6.17)  the  condensed  frequency 
function  of  the  zeros  of  g*(z)  —  0  is 

(8.13)  P{z) 

exp  {  —i  {q  -  Zy)H2Vft)~1(q  ~  Zy)  ]dg'  . 

Now  identify  the  matrix  Zvft  with  the  matrix  A  of  lemma  8.1,  taking  the  partitioning 
such  that  Au  is  the  variance-covariance  matrix  of  the  real  and  imaginary  parts  of  g*. 
Accordingly  we  take 

(8.14)  fli=  —  zoy  ,  (*2  =  ff'  —  ziy  , 

so  that  lemmas  8.1  and  8.2  yield 

(8.15)  P(z)  =  4^  Un  1 _,/2  \A«  -AuATMu  I  _,/7  («.+  *.T)  *  (o,+  iir) 

•  exp  {  —  ^aiAn1ai  -  i  (a2  —  A2iAn1ai)  ^  (Ajj  —  A2iAu1Ai2) 

•  (02  — A2iAn1ai)  }  da% . 

Put  b  =  a2  —  A2iAn1ai.  Then  (8.6)  gives,  using  results  from  Turnbull  and  Aitken  (see 
pp.  175-176  in  [17]), 

<8-i6>  p^=expLl 

•/exp  {  —  \b^B22b}  (6+ 2iY+A2iAn1ai)t  (&+ 2iY+A2iAn1ai) 

(  1  t  a  >  trace  { D } 

2-^-pw. 


a2  =  g  —  ziy 
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where 

(8.17) 

D  =  2?22*  +  ( Z\y  -^A-nAnaf)  *  ( Z\ y  +  A^\A\\a\) . 

Since 

(8.18) 

An  =  zoF 2o  >  A12  =  A| i  =  zoF z\  ,  A22 =  ZiV z\ 

we  obtain 

(8.19)  D  = 

ZiV z\  —  ZiF  zo  ( z0F  zo)  -1  z0F  zi 

+  I  2iY  —  ziVzl  ( 2oP2o)  1 2«r H  2i  Y  -  ZiV' 2o  ( UV zl)  1  zorl  * 

=  [2,-  (*iK&(M'rir1M[V+rr,Hli-  ( 2xK 2o)  ( 2<y 2o) -1 2o] f . 

We  have  thus  established 

Theorem  8.1.  If  the  coefficients  c*  —  { a*,  a*,  •  •  • ,  a*,  b*,  b*,  •  •  • ,  b*\  of  the  equation 
(8.20)  «*(*)  («,*  +  »4*)(*  +  iy)'- 0 

J=0 


are  distributed  normally  about  a  vector  mean  y  with  nondegenerate  variance-covariance  ma¬ 
trix  V,  the  condensed  frequency  function  of  the  roots  of  g*(z )  =0u 


(8.21)  P(z)  =P(x,y)  =H2_ 


_  1 


(zoy) 


t/® 


( z0F  2o)  1  ( 2or) 


2irn 


ZoV  zj 


1/2 


trace  (  £t) , 


where 


(8.22)  z0=  (1,  2,  z2,  •  •  •  ,  zn) ,  Z\—  (0,  1,  2  z,  •••  ,  »zn_1) , 

and 

(8.23)  <  =  zi  —  ( ziV2o)  (2oFzo)  _1zo  , 


and 

(8.24)  Af -F+yt1 


w  the  matrix  of  second  moments  of  c*  about  the  origin. 

In  deriving  this  result  we  have  assumed  that  V  is  nonsingular.  However,  in  view  of 
theorem  3.1,  theorem  8.1  remains  valid  provided  2oF2j  is  nonsingular. 


9.  Normally  distributed  real  coefficients 

The  case  of  normally  distributed  real  coefficients  arises  when  c*  is  distributed  normally 
about  a  mean  y  =  { cd,  0}  with  variance-covariance  matrix 


(9.1) 


We  shall  assume  that  the  distribution  is  nondegenerate  in  the  sense  that  V  is  nonsingular. 

As  we  have  just  remarked  in  section  8,  the  condensed  distribution  of  the  zeros  is  given 
by  theorem  8.1  provided  20FlJ  is  nonsingular.  Accordingly  we  shall  prove 
Theorem  9.1.  If  U  is  nonsingular  and  |y|  ^  8  >  0,  theorem  8.1  holds. 

Let 

(9.2)  x0=  (1,  *2,  ”•  >  *"), 


j  1,  2,  •  •  •  ,  w . 
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Set  X=  {xo,  JCi , •  *  - ,  xn} .  Since  |JT|  =  1,  XUX t  is  positive  definite.  By  the  matrix 
form  of  Taylor’s  theorem 

(9.3)  Zo  -  (1,  iy,  i2y2,  •  •  •  ,  inyn)X  . 


Hence 

(9.4) 


-y  iJT 

yoX, 


where  y0  +  iyi  =  (1,  iy,  Py2,  •  •  • ,  inyn)  and  y0,  yi  are  real.  Since  n  ^  1  and  |y|  ^  5  >  0, 
neither  y0  nor  yi  is  null,  and  there  is  an  (n  1)  X  (n  +  1)  nonsingular  matrix  Y  whose 
first  and  second  rows  are  y0  and  y\.  Thus 


(9.5) 


Zo 


V  O' 

/ 

*t 

2o  = 

6 

0 

V 

rt 

y  0 


ViXUX'  yl 


y*XUX'  yl 
VlXUX'  y\. 


is  a  minor  on  the  leading  diagonal  of  the  positive-definite  matrix  YXUX^Y^  and  is  there¬ 
fore  positive  definite  itself. 

There  remains  the  interesting  case  y  =  0.  We  shall  show  that  the  real  axis  supports  a 
line  density  of  the  condensed  distribution  of  the  zeros  of  g*.  We  write  in  place  of  (9.1) 

(9.6)  K  =  (^  tjJ),  0  <«-.'£  (lx)  >/» 

and  proceed  to  evaluate 

(9.7)  f(x)  =  lim  f  P (z) dy , 

v — >0+  J — $ 


which  will  then  represent  the  required  line  density  on  the  real  axis,  in  view  of  theorem  3.1. 
Define 


(9.8) 

Put 

(9.9) 

Then 

(9.10) 

(9.11) 
Thus 

(9.12) 
and 

(9.13) 

where 

(9.14) 

Hence 


I  jk  1 kj  XjXk 
Ui„=  Uhi  =  XjU  xl 

y  =  X  v* ,  |  y  |  ^  8  =  v1 ,  |  X  [  ^  v*1  . 


j,  k  0,1 ,  •  • '  ,  n  . 


ZoV 


20= ( 

*■  =  ( 

tl-( 


x0  —  X2  v16x2  , 

X  v*xi  , 

Xi  —  3X2  w16X3  , 


Xo 


2Xt>8x2  , 

Uoo-  2X2i>16£/02  , 
\v*Uoi 


Xi 


-2\v°x,  \+0(  21) 
-3\2v'«xJ  ^  K  } 

X  Uoi  \+Q(v2l\ 

|  20  P  Zo  |  =  A  d16  +  0(  a21) , 

A  =  tfoo/oo+  X2  (  UooUn  ~  Uh)  . 


(9. 1  5)  (  20Vzt)  -1  =  A-1  * 


\2U  n, 


—  X  v~8Uoi 


8  Uoi  ,  I-'16  Uoo  -  2  X2  f/02)  + ° (  V&) ' 
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It  is  evident  that  A  >  0  since  I70o  and 
(9.16)  U*1' 


(Un  U  oA 
\U01  Uj 


i  =  0,  1,  2,  •••  ,  n 


U oi  I/i 

are  minors  on  the  principal  diagonal  of  XUX t.  Next  write 

(9.17)  Aj-Xi a, 

From  (9.10) 

(9.18)  lot  =  (X°x»X»XA)+0(”!1)  ' 

Some  straightforward  algebra  now  gives 

(9.19)  -i(faT),(*,Vilr1(**r) 

=  (/oo^S  +  X!(C/„^-2C?o,^yl,+  (;oo^?))+0(tl1)> 

(9.20)  trace  (<Mc')  =  A_1Ioo ( Urn Un  -  Uh)  +A~2(  tWr  -  UnAo)*  +  <H  i>1) . 
Hence,  from  (8.21) 

/•noi\  r»/-x  1  i  loo  (  UooUn  ~ ’  Uh)  ,  lUUo0Al-UolA0)2} 

(9.21)  p(2)=^_5j - — - + - ^ { 

•exp  j-^[/oo^+X!(y„^S-2£/„1r4o^i+J/ooci;)]|+0(t-"’). 

The  exponential  term  in  (9.21)  is  never  greater  than  unity,  and  the  remaining  terms  in¬ 
volving  X  are  of  the  form  A“3/2  and  A~6/2.  We  have  to  integrate  over  —5  =  — 1>7  <  y  — 
Xt)8  <  »7  =  5.  Hence  we  may  multiply  (9.21)  by  v*  and  integrate  over  X  from  —  tr1  to 
tr1.  Allowing  v  — *  0  +  we  obtain  from  (9.7) 


(9-22>  i 


”  ( Zoo ( Um Un  —  Uoi)  ,  liar* A,-  U»A,y 


I 


A8/2  A6/2 

exp  j-^UooAl  +  \!(UnAS-2Uo,AoA1+UooA!)ll  d\. 


In  (9.22)  make  the  substitutions 
/oo£/oo 


(9.23)  M  = 
We  get 


X2 


«  TT  TT  ri2  „  ( UooAi  — ■  UoiAo) 2 

P=UooUn—Uoi,  q  — - y- - 


(924)  t(x)-ex vkz±AaU!™Lr  1  \t  i  qii  l 
(9.24)  /(*)  2wnUoo  Jo  ((*  +  #)«  («+*)) 

**• 


ft  *  Px/ 2  exp  (  -\A\!  Uqo) 

n)  =  irnUoo 


If  q  =  0,  (9.24)  gives 
(9.25) 

If  q  0,  we  make  the  substitutions 


(9.26) 
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and  find 

(9.27)  /W-^?lS,fr!lll  +  t(,)1- 
where 

(9.28)  K(r)  =  K/J  /  [*2-  (1  +  r2)]  e~*^ dt. 

rl  JT 

The  Hermitian  integrals  (see  p.  590  in  [11])  are  defined  by 

(9.29)  m-(f)  “/"* 

Thus 

(2  /ir")  l/2 

(9.30)  K{r)  =  -~/y7-  (2Hh2(r)  +2r  Hhi(r)  -Hh0(r)  } 


(2/t)  1/2 


Hhi  (r) , 


in  view  of  the  recurrence  relation 

(9.31)  (»+l)Hh,+1(r)  +rHhro(r)  -Hh,-i(r) 


0  . 


We  can  now  combine  (9.25)  and  (9.27)  into  a  single  result  and  collect  the  results  in 
Theorem  9.2.  If  the  coefficients  a*  =  (a*,  at,  *  •  • ,  a*)  of  the  equation 


(9.32) 


g*(z)  =2  a*z}'=  0 

J-0 


are  real  and  distributed  such  that  the  joint  distribution  of  g*(x)  and  dg*{x)/dx  is  a  non¬ 
degenerate  normal  distribution  for  each  fixed  x,  the  condensed  distribution  of  the  roots  of 
g*(z )  =  0  has  a  line-density  frequency  function  on  the  real  axis 


(9.33) 
where 

(9.34) 

(9.35) 


and 

(9.36) 


A0=EIg*(x)],  Ai  —E  [ , 
Uoo  =  var  [g*(s)]  ,  U01  -  cov  [g*(s) ,  * 


U 


11  =  var  [ 


dg*  (x)~| 

dx  J 


I  U00A1—  U01A0 1 
Ww(UwUn- Ul))^ 


Using  the  British  Association  tables  [1],  we  may  calculate  table  I.  When  r  is 
sufficiently  large  for  K(r )  to  be  negligibly  small  in  comparison  with  unity,  we  have 
from  (9.27) 


gl/8exp(-Mo/U00)  .  _  1 

(2t)  »/*(/«,  ax  V2x 


e  ri/2dT 


t  —  Aq  U  o/’ 


(9.37) 
where 

(9.38) 


nf  (x)  dx 


no 
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is  the  expected  value  of  g*  divided  by  its  standard  error.  It  will  then  be  a  reasonable 
approximation  to  take  the  percentile  points2  of  the  real  roots  of  g*(x)  as  the  real  roots 
of  the  corresponding  percentile  points  of  g*(x).  This  fact  was  previously  discovered  by 
Geary  [9]  for  the  particular  case  n  =  1.  Also  in  the  particular  case  n  =  1,  (9.37)  forms 
the  basis  of  Fieller’s  theorem  [6]  which  is  much  used  in  bio-assay.  In  the  usual  textbooks 
on  bio-assay  (for  example,  see  pp.  27-29  in  [8]),  there  seems  to  be  no  allowance  for  the 
fact  that  (9.37)  is  an  approximation  involving  the  neglect  of  K(r).  Fieller  and  Finney 
use  fiducial  theory.  I  do  not  know  on  what  logical  foundations  fiducial  theory  rests,  but 
I  can  imagine  that  Fieller’s  theorem  might  be  correct  within  a  framework  resting  on  such 
foundations  (if  they  exist).  The  particular  case  when  the  coefficients  a*  are  independ- 

TABLE  I 


f 

K(r) 

K(r ) 

0.0 

00 

1.2 

0.0935 

0408 

0.1 

7.0187 

0662 

1.4 

0.0523 

8306 

0.2 

3.0684 

4636 

1.6 

0.0290 

5246 

0.3 

1.7784 

0828 

1.8 

0.0158 

6176 

0.4 

1.1521 

9418 

2.0 

0.0084 

9070 

0.5 

0.7911 

8623 

2.5 

0.0016 

0331 

0.6 

0.5622 

4244 

3.0 

0.0002 

5477 

0.7 

0.4082 

2679 

3.5 

0.0000 

3342 

0.8 

0.3005 

1809 

4.0 

0.0000 

0357 

0.9 

0.2231 

8031 

4.5 

0.0000 

0031 

1.0 

0.1666 

3094 

5.0 

0.0000 

0002 

ently  and  normally  distributed  about  zero  mean  with  unit  variance  is  of  interest.  We  have 
U  —  I  and  a  =  0,  so  A0  =  Ax  =  0  and  (9.25)  gives 


(9.39)  f(x)  =1L°9*X_L.  * 

7 XtlL  00 

_  (1-fci) 

7T»(1  ~  X2)  ’ 


where  hn  is  defined  by  (1.  1).  Kac’s  result  [12]  follows  at  once  on  multiplication  by  n  and 
integration  over  —  <»  <£  x  5s  <» .  Actually  (9.39)  contains  more  information  than  Kac’s 
formula,  since  it  not  only  affords  the  average  number  of  real  roots  but  it  also  shows  how 
they  are  distributed  along  the  real  axis. 
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